Stationary and dynamical properties of a zero range process on scale-free networks 
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We study the condensation phenomenon in a zero range process on scale-free networks. We show that the 
stationary state property depends only on the degree distribution of underlying networks. The model displays a 
stationary state phase transition between a condensed phase and an uncondensed phase, and the phase diagram 
is obtained analytically. As for the dynamical property, we find that the relaxation dynamics depends on the 
global structure of underlying networks. The relaxation time follows the power law x ~ L z with the network 
size L in the condensed phase. The dynamic exponent z is found to take a different value depending on whether 
underlying networks have a tree structure or not. 
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I. INTRODUCTION 

A condensation is an intriguing phenomenon observed in 
nonequilibrium systems [1]. In interacting particle systems, 
a finite fraction of particles may be condensed onto a single 
site. Phase separation in driven diffusive systems 0001 or 
jamming in traffic flows 10] is also an example of the conden- 
sation. The condensation has been studied within the frame- 
work of a zero range process (ZRP) introduced by Spitzer |5]. 
Those studies revealed that the condensation can be caused 
by effective attractions between particles or quenched disor- 
der 00000 HolUlll . In the former the condensation is 
a spontaneous symmetry breaking transition. In the latter it 
is similar to the Bose-Einstein condensation occurring in the 
ideal Bose gas lfT3l . 

The condensation has been studied mostly on periodic lat- 
tices. On the other hand, recent studies show that lattices (net- 
works in general), upon which physical systems are defined, 
may have more complex structure 1131 M-4 1 1 51 llol . Partic- 
ularly many real-world networks are found to be scale-free. 
That is to say, they are characterized with the power-law de- 
gree distribution 



(i) 



where Pdeg(k) is the fraction of nodes (or sites) with the de- 
gree k (the degree of a node denotes the number of neighbors 
connected to it) and y is called the degree distribution expo- 
nent. This is contrasted to a peaked distribution for periodic 
regular lattices. 

The power-law distribution implies that the SF networks 
have an inhomogeneous structure. The structural inhomo- 
geneity is an important factor in understanding dynamical pro- 
cesses in SF networks d d |U E3 El- In random 
walks, for instance, a diffusing particle visits highly connected 
nodes more often than others. Moreover diffusion processes 
are asymmetric. It means that the speed of random walks from 
a node to another is not equal to that in the reversed direc- 
tion in general. The asymmetry can be quantified with the so- 
called random walk centrality 1 17], which is defined on each 
node and proportional to its degree. This example shows that 
each node plays a different role in stochastic processes in in- 
homogeneous networks. It also suggests that interactions in 



many particle systems would be important since particle den- 
sity could be high at some nodes. 

In this paper, we study a ZRP on SF networks and inves- 
tigate the interplay of the dynamics of interacting particle 
systems and the inhomogeneous structure of underlying net- 
works. The present work is an extension of our recent work 
presented in Ref. 12311 . which shows that the structural inho- 
mogeneity gives rise to the novel type of condensation in the 
ZRP. We will present a detailed review on the results. In ad- 
dition, we will elaborate a scaling theory for the relaxation 
dynamics in the condensed phase. Main results are following. 
The relaxation dynamics has a hierarchical nature; starting 
from small degree nodes, the relaxation proceeds to larger and 
larger degree nodes. The relaxation time x follows a power- 
law scaling x ~ Lr with the dynamic exponent z and the net- 
work size L in the condensed phase. The relaxation dynam- 
ics depends on the global structure of underlying networks so 
that the dynamic exponent takes a different value depending 
on whether the networks have a tree structure or not. Usually, 
the degree distribution is considered to be important in charac- 
terizing the dynamics of complex networks. Our result shows 
that the global structure of networks is an important ingredient 
in understanding the dynamics. 

This paper is organized as follows. In Sec.|II| we provide a 
general review on stationary and dynamical properties of the 
ZRP. Then in Sec. [H]] we introduce the ZRP on SF networks 
and present its stationary state phase diagram. In Sec. lIVI we 
investigate the relaxation dynamics. We summarize and con- 
clude the paper in Sec.lVl 



II. ZERO RANGE PROCESS 

Consider a system of N interacting particles on a network of 
L nodes. The total particle density is given by p = N/L. Each 
node i can be occupied by any integral number of particles (no 
exclusion). With n, ■ = 0, 1 , 2, • • • denoting the occupation num- 
ber at node i, a microscopic configuration is represented with 
n= (m,re 2 ,--- 

A ZRP is specified by jumping rates {/?,(«)} and hopping 
probabilities { ?}<-,•}. With those, at each node i with n, > 0, a 
particle jumps out at the rate /?;(«,) and then hops to a node j 
with the probability Tj<-i. Note that the jumping rate depends 
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only on the occupation number at departing nodes. This is an 
important feature of the ZRP, which makes the exact station- 
ary state probability distribution available 1 26]. Recently, gen- 
eralizations with multi-species particles [24] or with multiple- 
particle hoppings 12511 are being considered. 

A particle interaction can be incorporated into the jumping 
rate. If particles are noninteracting, they move independently, 
hence p(n) °< n . An attractive interaction can be imposed 
with the jumping rate function p(n) that grows sub-linearly or 
decreases with n, or a repulsive interaction with a super-linear 
jumping rate function. In the context of data packet transport, 
the jumping rate pi{n) corresponds to the number of packets 
that can be processed by a node i having n packets per unit 
time. That is, the jumping rate function models the transport 
capacity. 

It is well known that the stationary state property of the ZRP 
is exactly solvable For later uses, we list useful formula 

below. The stationary state probability distribution P s tat(n) is 
given by 



fstat(n) 



1 



Z(L,N) f = \ 



where 



m =n 



CO, 



„/ii \Pi{n') 



(2) 



(3) 



for n 7^ and /, (0) = 1, and Z(L,N) is a normalization con- 
stant. Here, {co,} is the solution of coupled linear equations 



(4) 



Note that these equations are equivalent to those for the 
stationary state probability distribution for a single random 
walker problem with the hopping probability {7}^,}. The 
left (right) hand side corresponds to the outgoing (incoming) 
probability flux. 

Introducing an auxiliary fugacity variable z fl], one finds 
that the mean occupation number m/ = («,) is given by 



where 



m = z^- In Fi(z) , 
oz 



n=0 



(5) 



(6) 



The series is defined for |z| < z c with z c being the radius of 
convergence. The fugacity should be determined from the 
self-consistent equation 



P = TE m 'W 



(7) 



A ZRP may display the condensation transition. We in- 
troduce two instances. A simplest ZRP displaying the con- 
densation might be the one with the jumping rate pi(n) = 



(1 — 8,,, ) + ocS;,, on a periodic regular lattice with nearest 
neighbor hoppings 1 1 ] . This model has a quenched disorder so 
that the impurity node z'o has a low jumping rate a < 1 . Apply- 
ing the formalism, one can show easily that the condensation 
occurs for p > p c = a/(l — a), where a macroscopic conden- 
sate forms at the impurity site with the occupation number 
m, = (p — p c )L, while the occupation number at the normal 
nodes is given by wi«,- = p c . 

The condensation also occurs in a system without quenched 
disorder. For instance, consider a ZRP on periodic regular lat- 
tices with the jumping rate pi(n) — p(n) = l +b/n. A straight- 
forward analysis shows that the system is in the condensed 
phase when b > 2 and p > p t = p c {b) CI 01- In this case, 
a macroscopic condensate forms at a node selected sponta- 
neously. 

As for the relaxation dynamics, an exact theory has not been 
known yet. However, qualitative features may be understood 
from a scaling argument El EH EH ■ To be specific, we 
consider the relaxation dynamic of the above example with 
the jumping rate p(n) = 1 +b/n in the condensed phase. Be- 
ing distributed uniformly initially, particles form small clus- 
ters at nodes due to a stochastic fluctuation. It is believed that 
a macroscopic condensate emerges by successive coarsening 
processes of the small clusters. Those clusters exchange par- 
ticles, merge into larger ones, and grow until there remains 
a single macroscopic condensate. This scaling picture pre- 
dicts the power-law growth of the cluster size in time and the 
power-law scaling of the relaxation time with the system size 
as x ~ U. It is known that the dynamic exponent z depends 
on the dimensionality of the underlying lattices and bias in 
particle hoppings 
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HI. ZRP ON SF NETWORKS 

A recent study reveals that structural inhomogeneity of SF 
networks leads to the so-called complete condensation l23ll . 
Being compared to the usual condensation in periodic regular 
lattices [1], the complete condensation has the features that (i) 
it occurs at any finite value of the particle density, i.e., p c = 
0, and that (ii) the whole fraction of particles are condensed 
onto a few high-degree nodes. We review these results in this 
section. 

Consider a ZRP on a SF network of L nodes. The network 
is assumed to have a power-law degree distribution Pdeg(k) ~ 
in the interval ko < k < k max . We are interested in sparse 
networks with finite mean degree k = J dk kPa eg (k). Hence, 
the degree distribution exponent is assumed to be larger than 
2,y>2. 

The node with the maximum degree k max is called the hub. 
While the lower bound ko is a constant, the maximum degree 
&max scales with the network size L. It can be estimated from 
the condition L P^ eg (k)dk — 1, which yields that fTsTl 



L p with p = 



1 



Y-l 



(8) 



Some SF networks having an explicit cutoff for £ max may have 
a different value of P [31]. In this work, however, we only 
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consider the generic SF networks with p given by Eq. 
For the jumping rate, we take 

Pi(n) =« 8 

for all nodes with a parameter 8 > 0. We take 

. — , if i and j are linked 

lj<-i = \ h 

, otherwise 



for the hopping probability (kj is the degree of a node z). That 
is to say, a particle out of i hops to one of its k( neighbors 
selected at random. 

First, one needs to solve the random walk problem in 
Eq. 10} for the stationary state probability P stat (n). The ran- 
dom walk problem with the hopping probability d 1 Oft was 
studied in Ref. 1 17], from which we find that 1 30] 

(Oi = k. (11) 

Inserting p(n) and {co,} into Eq. Q, one finds that 



\, Such a problem has been studied in the context of a disor- 

dered ZRP and a disordered exclusion process in one dimen- 
sion 01 0, [HI El- In those studies, the jumping rate distribu- 
(9) tion is taken to be Pj U mp(c < p' < 1) ~ (p' — c) v with the model 
parameter v and c for the minimum jumping rate. The present 
model at 8 = coincides with that with c = and v = y— 2. 

It is noteworthy that the mean occupation number in the 
stationary state is solely determined by the degree, which is 
(10) a local quantity, and independent of any other characteristics 
of a network. This, however, is not the case for the relaxation 
dynamics studied in Sec. lIVI 

At 8 = 1, the particle jumping rate is directly proportional 
to the occupation number. Particles move independently and 
the system reduces to a non-interacting system of N random 
walkers. So, the stationary state occupation number distribu- 
tion becomes proportional to the single walker distribution in 
Eq. il Q , mi CO, = kj. This solution is also obtained from the 
general formula in Eq. (114-1 . Since J\(x) = e*, m,-(z) = zki, 
which is proportional to the degree as expected. 

For general values of 8, the solution is non-trivial. We will 
solve the problem for 8 = and 8^0 separately. 



i=i 



which further yields from Eq. (|6ji that 

(zki)" 



F(z) = E 



as 



(12) 



(13) 



Hence, the mean occupation number in the stationary state is 
given from Eq. (0 by 



m,-(z) = x 



3 In Jg(x) 



dx 



where 



„to(»!) 8 



(14) 



(15) 



The fugacity z has to be determined from the self-consistent 
equation. 

We add a few remarks on the stationary state probability 
distribution. Although the jumping rate is the same at all 
nodes, the probability factor /;(n) in Eq. dl 2i has the node 
dependence due to the structural inhomogeneity of underly- 
ing networks. On the other hand, consider a ZRP in periodic 
lattices with disorde red jumping rates p\{n) = n 5 /ki and with 
nearest neighbor particle hoppings. It is evident that the two 
models share the same stationary state. That is, the structure 
inhomogeneity plays the same role as the quenched disorder 
in the jumping rate. At 8 = 0, in particular, the correspond- 
ing ZRP in periodic lattices is characterized by the disordered 
jumping rate p' = l/k whose distribution is given by 



dk 
dp 7 



(16) 



A. S = case 

At 8 = 0, the infinite series for Ts(x) has a closed form 
expression !Fo(x) = 1/(1 — x), which yields 



zki 
I -zki 



(17) 



The fugacity is in the range z < z c — 1 / ^max <C 1 ■ 

Note that m,-(z) is divergent for the hub as z approaches z c . 
So, it is convenient to decompose the self-consistent equation 
as p = p. s + p„ where 



Whub 



1 Zk n 



L LI Z&max 

is the particle density at the hub and 



1 zki 
u i^hub 1 



zki 



(18) 



(19) 



is the remaining density. For large L, one can evaluate p„ as 



l dk P ** m =zl- 1 
l—zk 



I'm 
*0Z 



dx 



(i 



The integrand becomes singular at both integral endpoints 
as z£ m ax — ► 1 and zko — > 0. Evaluating the singular parts, 
we obtain that p„ = o(z) + o(z Y_1 ln(l — zk max )), where the 
first (second) term is a contribution near the lower (upper) 
endpoint 13211 . Since z < l/k max — <C 1 and y > 2, p„ 
vanishes in the large L limit. Therefore we conclude that 
the whole fraction of particle should be condensed into the 
hub [complete condensation], 



P. S = »?hub/£ = P • 



(20) 
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The fugacity is then given by 

z=*«L/(i+i/(p£))' 



■L-P(l-l/(pL)). (21) 



So the mean occupation number at other nodes with <C £ n 
is given as 



iyhub — zki — ki/k n 



(22) 



which is proportional to the degree with the vanishing coeffi- 
cient. 



B. S > case 

For 8 > 0, the infinite series in Eq. J15i does not have a 
closed form expression except for Js=i i x ) = We approx- 
imate the series in the following way. Putting the summand 
as x"/(n!) 8 = exp[g(ra)] and using the Stirling formula, one 
obtains that 

g(n) ~ n\nx— 8 / — ln(27tn) +n\nn — n 
This function is peaked at n = no and behaves near n ~ no as 

g(«)~g(«o) + ^"(«o)(n-no) 2 + -" • 

The peak position is easily found to be no ^ x 1 / 5 from the 
condition g' (no) =0, which yields g"(no) = — 8/no + o(nQ 2 ). 
We then approximate the series with the Gaussian integral as 

J 5>0 (x) ~ [°° dn'exp[g"(n )n' 2 /2] 

J -no 

~ e 4 ' ( " o) /° O t/n'exp^"(n )n' 2 ] 



-^f2^V8 



(l-S)/2 



exp 



(23) 



In the second step, we extend the integration interval. This 
approximation is valid when |g"(no)|«o = S* 1 ^ 5 ^ U that is 
x>(l/8) 8 = o(l). For jc -c 1, one can simply approximate 
T§{x) with a few lowest order terms as Ts(x) = 1+x+O (x 2 ). 

Note that the approximation in Eq. (I23i coincides with the 
exact result at 8 = 1 . We have tested its validity at other values 
of 8 numerically. Figure ^ shows that the approximation is 
very good unless x ~ 0. 

Using the approximation for T§(x), the mean occupation 
number is given by 



zki for zki < 1 

(zki) 1/s for zki>\. 



(24) 



The mean occupation number of a node increases monoton- 
ically with its degree. So we expect that nihub ^ 1 an d that 
zkmax. > 1 ■ On the other hand, for a node with the minimum 
degree Icq, one may have zh) <C 1 or zko > 1 depending on the 
magnitude of z- Hence, one must consider the following two 
cases separately. 




FIG. 1: Comparison of f$(x) (circular symbols) defined in Eq. I15i 
and the approximation (solid line) in Eq. <23t at 8 = 0.5 (a) and 
1.5(b). 



(i) Firstly, let us assume that the fugacity is in such a range 
that zki > 1 for all nodes. This is the case when z is a non- 
vanishing constant in the large L limit. Then the occupation 
number is given by m, ~ (fc/z) 1 / 8 for all nodes, and the self- 

consistent equation has the solution z — p 5 / |V/ 8 ^ where 
/t 1 / 8 = // o max dk k l / & P deg (k). It could be a non-vanishing con- 
stant only when A: 1 / 8 ~ fb!** dkk~^~ 1 /^ remains finite in the 
large L limit. It imposes that 8 > 8 C with 



y-1 



In this regime, we find that 



rrij ~ ck. 



1/8 



(25) 



(26) 



for all nodes with a constant c = p/ [k l ^j. The occupa- 
tion number at the hub with the degree k max ~ LP scales as 
m hu b ~ i p/S - Since p = 1 /(y- 1) and 8 > 8 C = 1 /(y- 1), the 
occupation number at the hub scales sub-linearly in L. There- 
fore, the condensation does not occur when 8 > 8 C . 

(ii) Secondly, let us assume that the fugacity z is in such 
a range that the crossover degree defined as k c = l/z scales 
with the network size L and is in the interval k$ <§; k c <§C k max . 
Then, the self-consistent equation becomes p = p„ + p s where 



(27) 



p„ =k c 1 / dkkP deg (k) 

Jko 



is the density at small-degree nodes with k < k c and 

C ^max 



P,=k^ 



' dkk x l h P Atg {k) 



(28) 



is the density at large-degree nodes with k > k c . The inte- 
gral part in Eq. |27i is smaller than the mean degree k that 
is finite (y > 2). So, p„ vanishes as ~ k^ 1 in the large L 
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limit, which yields that p s = p. In order to have a finite value 
of p s , J^* dkk l l^~7 should be divergent, which yields that 
8 < S c = l/(y — 1). And it should be of the same order as 
k c which yields 

f [lnfc ma x] 8c , for 8 = S C 
k c ~ <^ (29) 
I [fcmax] 1 - 8 ^ , for 8<S C . 

With the crossover degree k c = l/z, the occupation number 
distribution is given by Eq. J24> .which can be summarized 
with the scaling form 

m i = G 8 (k i /k c ) (30) 

with the scaling function Gg (y ) that behaves as Gg (y <C 1) ~ y 
and G s (y > l)-// 8 

We add a few remarks on the result in Eq. (13 0> . Nodes 
with ki < k c covers the whole fraction of nodes in a network 

since f£ max dkPd eg (k) ~ k c ^ ^ vanishes in the L — > °° limit. 
However, the total density of particles on those nodes, given 
by p„ ~ k~ , vanishes in the thermodynamic limit. There- 
fore, the whole fraction of particles are condensed onto the 
vanishingly small fraction of nodes with kj > k c . For those 
nodes, the occupation number is given by Wj ~ (ki/k c ) 1 ^ ~ 
(^maxAc) 1 ^ 8 (^i7^max) 1//5 - It nas the overall L dependence as 
(&maxA t ) 1/8 ~ i p/8r - L for 8 < 8 C . and - L/\nL for 8 = 8 C . 
That is, the occupation number scales linearly with the system 
size L for 8 < 8 C . One can also find that the scaling form in 
Eq. d30t can be applied to the case of 8 = studied in the pre- 
vious subsection with k c = £ max - Therefore, we conclude that 
there is the complete condensation for 8 < 8 C . At the threshold 
8 = 8 C , the logarithmic correction appears. 

As explained earlier, the stationary state of the present 
model is equivalent to that of the disordered ZRP on periodic 
lattices with the jumping rate p\(n) = n^/kj. The disordered 
ZRP with n-independent jumping rates, which corresponds to 
the 8 = case, has been studied in Refs. B8l l9l flQl . Those 
studies show that the condensation is triggered by the node 
with the slowest hopping rate and the complete condensation 
occurs when the slowest jumping rate vanishes IHIKa IKiII ill . 
Generalizing the results, one can understand our phase dia- 
gram intuitively. The slowest hopping rate node corresponds 
to the hub with the maximum degree k max ~ zA Assuming 
that the node is occupied macroscopically, one finds that the 
jumping rate is given by p' = n 8 '/ki ~ L S ~P. It vanishes when 
S<P=l/(8— 1), which coincides with the phase boundary. 

Summarizing the results, we present the stationary state 
phase diagram of the ZRP on SF networks in Fig. [2] The 
region I corresponds to the uncondensed phase, where the oc- 
cupation number is given by Eq. J26I . The region II corre- 
sponds to the condensed phase, where the occupation number 
is given by Eq. (I30i . The phase boundary is given by 8 = 8 C 
with Eq. d25i . The phase diagram suggests that the condensa- 
tion also occurs in random networks (y = °°) at 8 = 0, which 
is observed in a numerical study [ 33] . The occupation number 
distribution in each phase is sketched in Fig.|2] 




FIG. 2: Stationary state phase diagram. The region I and II corre- 
sponds to the uncondensed and condensed phase, respectively. 



IV. RELAXATION DYNAMICS 

The preceding analysis shows that the occupation number 
distribution in the stationary state is solely determined by the 
degree distribution Pdeg(k) being independent of any other 
characteristics of underlying networks. However, the relax- 
ation dynamics may be affected by the structure of underlying 
networks. We are interested in the relaxation dynamics in the 
condensation phase. Unlike the stationary state property, an 
exact theory is not available for the dynamical property of the 
ZRP. We rely on Monte Carlo simulations for the study and 
provide a scaling argument to understand dynamical scaling 
properties. 

We describe our Monte Carlo simulation algorithm. Ini- 
tially at time t =0, N particles are distributed uniformly and 
randomly on a network of L nodes. At each step, a parti- 
cle is selected at random. Let i be its residing site and n; 
be the occupation number there. Then, with the probability 
p(rij) /rii = nf _1 , the particle jumps out of i and hops to one of 
its neighboring nodes. With the probability 1 — the 
selected particle does not move. The time is measured in unit 
of one Monte Carlo sweep consisting of N trials. Note that the 
particle-based dynamics with the modified jumping probabil- 
ity p{n)/n is equivalent to the original node-based dynamics 
with the jumping rate p(n). 

As for SF networks, we use the Dorogovtsev-Mendes- 
Samukhin model 1 34] that is a generalization of the Barabasi- 
Albert model 1 35]. With the model, one can generate networks 
of tree structure or looped structure with arbitrary values of the 
degree distribution exponent y. 

In order to get an insight into the nature of the relaxation 
dynamics, we have observed the time evolution of the mean 
occupation number distribution. The Monte Carlo simulation 
data are presented in Fig. [3] (a). They are obtained on the 
looped SF network of L = 10 4 nodes with y = 3 (8 C = 1 /2). 
The particle density is taken to be p = 1. We take the jump- 
ing rate p(n) — n 5 with 8 = 0.2 so that the system is in the 
condensation phase. The plotted quantity is the ensemble 
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FIG. 3: (a) The occupation number distribution at < / < 2 15 . 
The dotted (dashed) line has the slope 1 (1/8 = 5). (b) Schematic 
plot of mft vs k in the log-log scale in the transient time. 



average of the occupation number n, averaged over nodes with 
the degree ki = k. 

The data show how the distribution, starting from the uni- 
form distribution at t = 0, approaches the stationary distribu- 
tion. One can see that the initial uniform distribution evolves 
very quickly (at t = 1) into a linear distribution ~ k. This 
is due to the random walk nature of the initial dynamics. The 
initial uniform distribution implies that the jumping proba- 
bility is the same for all particles at t = 0. So they behave 
as non-interacting random walkers, which has the distribu- 
tion rrik ~ k. After that, the distribution evolves further and 
approaches the stationary distribution. The MC data shows 
that the system has already reached the stationary distribution 
given in Eq. J30l > at t = 2 15 . 

In the transient period, we find that there are three differ- 
ent regimes that are separated by two time-dependent degree 
scales denoted by ko and %. It turns out that the occupation 
number distribution behaves in each regime as 




for k<kp 

for ko < k < ku 

for kjj <k . 



(31) 



We illustrate the behavior with the schematic plot in Fig. 13(b). 
The degree scales grow in time and approach ko — » k c with 
the crossover degree k c given in Eq. j29t and ku — > k max in the 
stationary state limit. Furthermore, the behavior of m<- for k < 
ku is identical to the stationary state distribution in Eq. d30i 
with ko and % playing the role of k c and £ max , respectively. 
For nodes with k > ku, °< k. They are still in the random 
walk regime. 

These observations lead us to conjecture that the relaxation 
dynamics has a hierarchical nature: A sub-network of small- 
degree nodes with k < %(?) is equilibrated first and large- 
degree nodes become equilibrated as the degree scale ku(t) 
grows in time. In a transient time t, those nodes with k < 
ku(t) reach the stationary state of a smaller network of size 
U(t) ~ with p = l/(y- 1). The other degree scale 

ko(t ) plays the role of the crossover degree scale in the smaller 




FIG. 4: Illustration the hierarchical relaxation dynamics. Up to an 
intermediate time step t, the subnetwork of L'(t) nodes and with the 
largest degree ku(t) ~ L'(t)$ reaches the stationary state. The equi- 
librated area expands as ku increases until it reaches k m - dx . 



network, that is, 



■l-5/5 c 



(ln%) ? 



for 8 < 8 e 
for 8 = 8 C 



(32) 



following from Eq. i29\ . This hierarchical nature is illustrated 
in Fig.|U 

The conjecture predicts that the relaxation dynamics is de- 
termined by the temporal evolution of the degree scale %(f ). 
We develop a scaling theory for ku(t). Suppose that the sys- 
tem is equilibrated up to a sub-network of L' nodes with the 
maximum degree ku ~ L'$ with p = l/(y— 1). The node with 
the maximum degree will be called a temporary hub. In or- 
der to proceed to the next hierarchy, particles should be trans- 
ferred to higher degree nodes. We will assume that particle 
transfers from the temporary hub to the next hub dominate the 
relaxation time scale, which will be confirmed numerically 
shortly. Since the sub-network is in the stationary state, there 
are n ~ L' ~ k\j particles at the temporary hub. During the 
unit time interval, An — p(n) ~ L' 5 particles jump out of the 
temporary hub and perform the random walk motions. They 
may wander within the sub-network and return to the tem- 
porary hub, or escape to the higher degree node in the next 
hierarchy level. Let P esc (L') be the escape probability for the 
latter process. Then, we obtain the following scaling relation 



1 / n 

(Z7) x VA^ 



xL 



n-5 



(33) 



The escape probability depends on the structure of the un- 
derlying networks. Suppose that the underlying SF network 
has a tree structure with no loop. On a tree a path from one 
node to another is unique. So, particles can escape from the 
temporary hub only when they hop along the link to the right 
direction toward the next hub among % links. Hence the es- 
cape probability is given by 



k- 1 
K U 



r/-p 



(34) 



for a tree network. On the other hand, in general networks 
with many loops, there are multiple pathways between two 
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nodes. Hence, we expect that the escape probability does not 
scale algebraically with L' or %, that is, 

Pe S c(i/)~i ,0 ~4 (35) 

for a looped network. Combining Eqs. (I33i d34i . and J35l >. we 
obtain that 

L ^ f l/(l+p-8) or % ^ f P/(l+P-8) (36) 

for a tree network and 

LWa-S) or kv-t^-V (37) 

for a looped network. The scaling behavior of Icq is obtained 
from Eq. d32i . 

The system reaches the global stationary state when L' be- 
comes equal to the network size L or, equivalently, k\j be- 
comes equal to k max ~ zA The power-law scaling of L' and 
£{/ implies that the relaxation time also scales algebraically as 

x ~ L z . (38) 

The dynamic exponent z depends on the global structure of 
underlying networks and is given by 

_ J 1 + p — 8 , for a tree network 

1 1 — 8 , for a looped network. " ' 

At the transition 8 = 8 C -, there is a logarithmic correction to 
the power-law scaling. 

We add a remark on the result in Eq. (1391 . As explained in 
Sec. m] our model at 8 = has the same stationary state as 
the disordered ZRP in periodic lattices with the jumping rate 
distribution Pjumpip') = p' v with v = y — 2. The relaxation 
dynamics of the ZRP has also been studied in Refs. OUGil]. It 
is known to have z = (v + 2)/(v+l) with the biased hopping 
in one dimension lattices. Interestingly, it coincides with the 
result for the ZRP on tree structure SF networks (p = 1/ (y— 
1 )). In the disordered ZRP the relaxation proceeds via a lateral 
growth of a domain within which particles are condensed on 
to the node with the slowest jumping rate. This mechanism is 
identical to the hierarchical relaxation in the SF tree networks. 

We verify numerically the analytic results. We have per- 
formed the Monte Carlo simulations on SF networks with var- 
ious values of y and with several values of 8. In order to 
study the structure dependent relaxation dynamics, we have 
performed the simulations both on tree networks and looped 
networks. 

According to the scaling argument, the degree scale %(f) 
separating the second and the third region in Eq. ( I32l > scales as 
ky[t) ~ t P/ z . And the occupation number at the node with k = 

kjj scales as ~ k\j^ ~ t l l z . Hence, the existence of % can 
be verified from the scaling plot of mk(t)/t^l z against k/t$l z 
at transient times t <C X. In Fig. 0(a), we present the data 
for the ZRP with 8 = 0.2 on a tree-structured SF network of 
L = 10 5 nodes and with y = 4. All data at different t collapse 
very well near k ~ $l z with the proposed value of the dynamic 
exponent z = 1 + P — 8. Figure|5](b) shows the data of the ZRP 




-2 2 -202 

10 10 10" 10" 10 10" 




FIG. 5: (color online) Data collapse analysis of Monte Carlo data of 
the ZRP with 8 = 0.2 on SF networks of L = 10 5 nodes and with y = 
4. Data in (a) and (c) are from the SF networks with the tree structure, 
while those in (b) and (d) are from the looped SF networks. The data 
are measured in the transient period at / = 4 2 (black), 4 3 (red), • • • , 
4 6 (magenta). 

with the same value of 8 = 0.2 on a looped SF network with 
the same value of y = 4 (8 C = P = 1/3). We also observe the 
perfect data collapse near k ~ $l z with the proposed value of 
the dynamic exponent z = 1 — 8. These support the structure 
dependent dynamic scaling of the %. 

Our scaling argument also predicts the existence of the 

degree scale ko(t) ~ k\j 8 ^ 5c which plays the role of the 
crossover degree k c for nodes with k < kjj(t). The scal- 
ing ku(t) ~ ftl z implies that k D ~ f(P- s PA0/z. Note that 
P = S c =l/(y — 1). Hence, it can be verified with the scaling 
plot of m.k against k/t^~^l z . Since the occupation number 
at the crossover degree is a constant, the occupation number 
need not be rescaled. We present the scaling plots in Fig.[5](c) 
and (d) for a tree SF network and a looped SF network with 
the corresponding dynamic exponent, respectively. Data from 
different times t collapse very well near k ~ t (P~ 5 V Z , which 
supports the analytic result. We have also performed the data 
collapse analysis at several other values of y and 8, and ob- 
tained the consistent result with the analytic result. 

In addition to the data collapse analysis, we also measure 
the dynamic exponent z through a temporal scaling behavior 
of /«hub- The scaling picture predicts that nodes with k < % 
are equilibrated while nodes with k > kjj are in the random 
walk regime with « k. So the occupation number at the 

hub scales as ~ k\j^ x (k max /ku). Since % ~ fil z in the 
transient times t -C x, the occupation number at the hub shows 
the temporal scaling 

m hu ^t^ Uz . (40) 

From the temporal scaling of mi lu b, we obtained the values 
of the dynamic exponent z at different values of y and 8 on SF 
networks of size L — 10 5 of the tree and looped structure. We 
present the result in Fig. [6] All the results are consistent with 
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FIG. 6: Dynamic exponent of tree (filled symbols) structure SF net- 
works and looped (open symbols) structure SF networks. Each sym- 
bol stands for y = 5 /2 (circle), y = 3 (square), and y = 4 (diamond). 
The solid and dashed lines represent the analytic result in Eq. 1391 
for tree and looped networks, respectively. 

the analytic results in Eq. (I39i . Approaching 8 = 8 C = 1/ (y— 
1), the numerical results seem to deviate from the analytic 
predictions. It is attributed to strong finite size effects near 
8 = 8 C for we obtained better results with larger system sizes. 

The numerical analysis supports the conjecture of the hier- 
archical relaxation and the scaling argument for the structure- 
dependent dynamic exponent z in the condensed phase with 
8 < 8 C . Our result shows that the dynamical property depends 
on the global structure of the underlying networks, although 
the stationary state property does not. 

V. SUMMARY 

We have studied the stationary and dynamical properties 
of the ZRP with the jumping rate p(n) = n s on SF networks 



with the degree distribution exponent y. The stationary state 
property is determined by the degree distribution of underly- 
ing networks. The model displays the stationary state conden- 
sation transition and the phase diagram is given in Fig. [2] The 
phase transition is driven by the structural disorder in underly- 
ing scale-free networks. It is distinct from the Bose-Einstein 
type condensation transition found in systems with disordered 
jumping rates, and from the spontaneous symmetry breaking 
transition found in systems on periodic lattices. In the con- 
densed phase 8 < 8 C = l/(y— 1 ), the mean occupation number 
distribution follows the scaling form of Eq. d30i . It shows that 
the whole fraction of particles are concentrated on higher de- 
gree nodes with k > k c given in Eq. ( I29l l. On the other hand, 
the dynamical property depends also on the global structure 
of underlying networks. The relaxation dynamics is found to 
have the hierarchical nature that it proceeds from lower degree 
nodes toward higher degree nodes. It is characterized by the 
power-law growth of the degree scale kij. From the scaling 
argument, we show that it grows in time with the different ex- 
ponents depending on whether underlying networks has a tree 
or looped structure (see Eqs. fl36i and i3H). Consequently, 
the dynamic exponent for the relaxation time x ~ U also de- 
pends on the global structure of networks (see Eq. ( I39» . Our 
result shows that the global network structure is the important 
factor in understanding the dynamical property of complex 
networks. 
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